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A new constraint of clustering for the AMD calculation is proposed. This new constraint 
gives us large improvement in studying the cluster structure by AMD which sometimes meets 
difficulty in giving rise to some specific cluster configurations. The usefulness of this new 
constraint is verified by applying it to the the study of the third K* — + band of 20 Ne 
which has been discussed to have 2a- 12 C structure. This band has not been easy even to 
construct by AMD. We see that the AMD+GCM calculation by the use of the new constraint 
gives rise to the third K* = + band which contains the 2a- 12 C structure as an important 
component. 

§1. Introduction 

For the study of the coexistence of the cluster structure and the mean-field-type 
structure in nuclei, the AMD (antisymmetrized molecular dynamics) method has 
proved to present us with a powerful approach for both stable and unstable nu- 
clei. 1 )' 2 ) AMD is a kind of ah initio theory which can describe any kinds of nuclear 
structure but does not require any model assumptions on the nuclear structure such 
as the assumption of the formation of the mean-field-type structure and the assump- 
tion of the existence of any clusters. One of the frequently used processes of AMD 
study of nuclear structure is the constrained AMD calculation in which we calculate 
the parity-projected energy surface as a function of some constraint parameter like 
in the constrained Hartree-Fock calculation. As the constraint parameter we usually 
use the quadrupole deformation parameter /3. If the system prefers to have a cluster 
state, it usually appears as the energy minimum point on the parity-projected energy 
surface as a function of the j3 parameter, which we abbreviate as the /3-constraint 
surface hereafter. However, the use of the /3-constraint surface can not be always a 
good way for the study of the cluster structure. It is because the quadrupole de- 
formation parameter is not a direct quantity to characterize the clustering structure 
but is a good quantity to characterize the mean- field- type structure. The minimum 
energy state for a given value of /3 tends to have more mean-field-type character than 
cluster structure character. Therefore the states on the /3-constraint surface tend to 
underestimate the clustering character. The state having more prominent clustering 
character is often located higher above the state on the /3-constraint surface at each 
value of /3. 

The inappropriateness of the constraint parameter /3 is sometimes rather serious 
because some kind of cluster structure which is reported to exist in experiments 
can not be easily obtained on the /3-constraint surface. One such example is the 
K w = 0™ band which has a- 40 Ca structure. 3 )' 4 ) As is reported in Ref.5), this K w = 
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CP band with a- Ca structure is not easy to obtain on the /3-constraint surface of 
AMD and hence in Ref.5) it was constructed by explicitly using the a- 40 Ca Brink 
wave function as the initial state for the frictional cooling process of AMD which 
is an energy variation process in AMD analogous to the imaginary time method 
in the Hartree-Fock approach. Another example is the third K n = + band in 
20 Ne which is strongly suggested by experiments to have 2a- 12 C structure since the 
band member states are strongly populated by the 8 Be transfer reaction on 12 C. 6 ) 
The assignment of the 2a- 12 C structure to this third K w = + band was supported 
by the semi-microscopic calculation of the ( 12 C- 8 Be) + (a- 16 0) coupled channel 
OCM (orthogonality condition model). 7 )' 8 ) However, as is seen in Ref.2), this third 
K w = + band of 20 Ne could not be shown even to exist on the /^-constraint surface 
of AMD. Hartree-Fock calculations have also never shown clearly the existence of 
the third K n = 0+ band of 20 Ne. 

It is clear that the above-mentioned deficiency of the /3-constraint surface for 
the proper treatment of the clustering problem is due to the inappropriateness of 
the constraint parameter for the study of the clustering. We should invent a new 
constraint parameter which is more suitable for treating the clustering character 
and should calculate the parity-projected energy surface as a function of the new 
constraint parameter. The purpose of this paper is to propose a new constraint for 
AMD which is suitable for treating the clustering problem and to verify its usefulness 
by applying it to the study of the third K w = + band of 20 Ne. We will see that 
we succeed in reproducing the existence of the third K 77 = + band of 20 Ne. We 
will also see that this calculated band contains the 2a- 12 C structure as an important 
component. 

In these days, the study of the coexistence of the cluster structure and the 
mean-field-type structure has become a rather hot issue for nuclei around the mass 
number 40. In this region of nuclei, recent experimental studies have assigned many 
superdeformed rotational bands which have the structure of many-particle many- 
hole configuration. 9 )~ n ) Many of the low spin states of these superdeformed bands 
have long been known to exist 12 ) and their natures have been understood often 
in relation with the clustering 13 ^ or quarteting 14 ) correlations. Therefore we are 
now confronting the situation which requires us to clarify the relation between the 
mean-field-type structure with superdeformation and the cluster structure with large 
deformation. We expect that our new constraint will play a powerful role in the study 
of the coexistence of the cluster structure and the mean-field-type structure in nuclei 
around the mass number 40. 

The content of this paper is as follows. In the next section, ^2 we explain 
briefly the formulation of the deformed-basis AMD 1 )' 2 )' 15 ) and in £j3]we introduce 
and explain a new constraint for clustering which we call the d-constraint. In £jU we 
prepare the basis wave functions to be used in our present AMD+GCM (generator 
coordinate method) calculation for 20 Ne. The basis wave functions consist of those 
constructed with the d-constraint in addition to those constructed with the usual 
/3-constraint. Two types of cluster configurations are adopted in constructing the 
basis wave functions with the d-constraint; one is the a- 16 type and the other is 
the 2a- 12 C type. In ^ 

we discuss the results of the GCM calculation. We will see 
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that our GCM calculation gives rise to the third K n = + band which contains the 
2a- 12 C structure as an important component. Finally in ^Hl we gi ye summary. 

§2. Framework of Deformed-basis AMD 

2.1. Wave Function and Hamiltonian 

In deformed-basis AMD, the intrinsic wave function of the system with mass A 
is given by a Slater determinant of single-particle wave packets, 

$ int = —rrf^Wu^, ■■■ , <Pa}, (2-1) 
VA\ 

Pi = (t>i( r )Xi£i, (2-2) 

where ipi is the ith single-particle wave packet consisting of spatial fa, spin Xi an d 
isospin £j parts. Deformed-basis AMD employs the triaxially deformed Gaussian 
centered at Zi as the spatial part of the single-particle wave packet: 



fa(r) oc exp 



E ■ 

'=x,y,z 



Zicr 



Xi = ct%X\ + PiXb Wi? + IAI 2 = !, 

£i = proton or nertron. (2*3) 

Here, the complex number parameter Z\ which represents the centroids of the Gaus- 
sian in phase space takes an independent value for each nucleon. The width param- 
eters v x , v y and v z are real number parameters and take independent values for each 
direction, but are common to all nucleons. Spin part Xi 1S parametrized by and 0i 
and isospin part £j is fixed to up(proton) or down(neutron). Z^, v x , z/y, v z and Oj, Pi 
are the variational parameters and are optimized by the method of frictional cooling 
explained in the next subsubsection. As the variational wave function, we employ 
the parity projected wave function: 

$± = P±$ int = i^$ int , (2-4) 

here P x is the parity operator and $i n t is the intrinsic wave function given in Eq. (|2-1|) . 
The Hamiltonian used in this study is as follows; 

H = f + V N + V c -f G , (2-5) 

where T and Tq are the kinetic energy and the energy of the center of mass motion, 
respectively. We have used the Gogny D1S force as an effective nuclear force V^. 
Coulomb force Vq is approximated by the sum of seven Gaussians. 

2.2. Energy Variation, Angular Momentum Projection and Generator Coordinate 
Method 

We perform the variational calculation and optimize the variational parameters 
included in the trial wave function Ea. l|2-4|1 to find the state that minimizes the 
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energy of the system E 



^ ($ ± |?i:|$ ± ) ~ 

2 = ,^+1^+, , H = H + V cnst . 



E± = ( g4|g±) \ n = H + V cnst . (2-6) 



We add the constraint potential V^nst to the Hamiltonian H to obtain the minimum 
energy state under the optional constraint condition. In this study, we constrain 
matter quadrupole deformation and impose a new constraint on the distance between 
quasi-clusters, by employing the potential V cns t = ^f n st((/^ 2 ) — /^o) 2 anc ^ u cnst((^ 2 ) ~~ 
cZq) 2 and we obtain the optimized wave function $^(/3o) = P ± ^ nt (Po) and <J>^(do) = 
P ± ^f Ilt (d ). Asforthedist ance between quasi-clusters, we will explain it in the next 
subsection. The evaluation of the quadrupole deformation parameter (5 is explained 
in Ref.16). At the end of the variational calculation, the expectation value of V^ ns t 
should be zero in principle, and in the case of /3-constraint we confirm that it is less 
than 100 eV. However, in the case of the constraint on the distance between quasi- 
clusters, the resultant (<i 2 ) value sometimes differs slightly from the aimed value 

(g?o), and hence the expectation value of V^ ns t, (v c nst^, becomes non-negligible. In 

such a case, we of course subtract (v cns t^j from E^ of Eq. ([2-6|) in calculating the 
energy of the system. 

Energy variation with the AMD wave function is performed by the frictional 
cooling method. The reader is referred to Refs.17), 18) for a more detailed expla- 
nation. The time development equation for the complex number parameters Z^on 
and Pi is as follows; 

dXi ndE ± . ,s , _s 

-lT = TidXf> (*=1,2,-,A) (2-7) 

and for the real number parameters v x , v y and v z ; 

dv a (jf dE ± 



dt h dv a 



(a = x,y,z) (2- 



Here Xi is Zi,ct{ or Pi. fi and // are arbitrary negative real numbers. It is easy to 
show that the energy of the system decreases as time develops, and after sufficient 
time steps we obtain the minimum energy state. 

From the optimized wave function, we project out the eigenstate of the total 
angular momentum J, 

*d J M K ( d o)=PMK*i(d ). (2-10) 

Here P(^ K is the total angular momentum projector. The integrals over the three 
Euler angles included in the P'lix are evaluated by numerical integration. 

Furthermore, we superpose the wave functions Q$ J MK and $dMK which have 
the same parity and angular momentum but have different values of deformation 
parameters Pq and distance parameters do and K. Thus the final wave function of 
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the system becomes as follows: 

^ = c n ^ K (f3 ) + c f n ^ K ,(/3' ) + ■■■ 

+d n <S> d J ± K „{d ) + d' n <f> d J ± K „,(d ) + • • • . ( 2 - n ) 

where quantum numbers other than total angular momentum and parity are rep- 
resented by n. The coefficients c n ,c' n ,--- ,d n ,d' n ,--- are determined by the Hill- 
Wheeler equation, 

8 ((^±\H\$i±) - e n ($ J n ±\$ J n ±)) = 0. (2-12) 

§3. New Constraint (d-constraint) 

In this paper, we propose a new constraint for creating the wave functions of 
cluster structure in AMD. This constraint is imposed on the distance between the 
centers of mass of quasi-clusters. The meaning of "quasi-cluster" will be soon ex- 
plained below. We call this constraint d- constraint. 

At first, we decide proton number and neutron number of each quasi-cluster. 
Next, we make a numbering of nucleons to fix which nucleon belongs to which quasi- 
cluster. By the "quasi-cluster" , we simply mean a set of nucleons which consists of 
given number of neutrons and protons. We do not impose any other properties to 
the quasi-cluster such as the spatial distributions of nucleons of the quasi-cluster. In 
this paper, we have studied the a- 16 and 2a- 12 C structures of 20 Ne. Therefore we 
treat quasi-a cluster and quasi- 12 C cluster and quasi- 16 cluster in 20 Ne. 

The center of mass of quasi-cluster C n , R n , is defined in the following way: 

Rna = -T- Re (a = x,y,z), (3-1) 
An iec n v" ff 

where A a is the mass number of quasi-cluster C n , and i £ C n means ith nucleon is 
contained in the quasi-cluster C n . It is to be noted that the a(= x, y, z) component of 
the spatial center of one particle wave function, is Rc ?i? . We define the distance 
d between quasi-clusters C n and C m , as 

(d 2 ) = \R n - R m \ 2 . (3-2) 

If the system favors to have clusters C n and C m with the mutual distance d = do, 
we will obtain, after the frictional cooling, the wave function in which A n nucle- 
ons gather closely around R n and A m nucleons gather closely around R m where 
\R n — R m \ = do, respectively. On the contrary, if the system does not favor to have 
clusters C n and C m with d = do, it can happen that A n nucleons and A m nucleons 
do not gather closely around R n and R m , respectively, where \R n — R m \ = do- Fur- 
thermore, even when A n and A m nucleons gather closely around R n and R m with 
|i? n — Rm\ = do, respectively, these A n and A m nucleons do not necessarily form 
the ground state configurations of the clusters C n and C m , but form in general the 
distorted configurations from the ground state configurations. In summary, the use 
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of ri-constraint does not necessarily mean to result in the construction of the C n -C m 
cluster structure, and even when the C n -C m structure is formed, the wave function 
created by ri-constraint is not necessarily the same as the Brink model wave function 
as usual cluster models. In the usual Brink model, nucleon wave packets of a cluster 
is at the same place, which means that the usual Brink model can only describe the 
cluster at approximate ground state. 

As examples of calculating 
2-body and 3-body cluster wave 
functions, we discuss "a- 16 0" 
and "2a- 12 C" wave functions 
which we use later in this pa- 
per (Fig^). The central and 
right figures are those calcu- 
lated by d-constraint. In get- 
ting the "a- 16 0" wave func- 
tion of the central figure, the 
constraint distance d between 
quasi-a and 16 O clusters is set 
to d = 4.5 fm, while in getting 
the "2a- 12 C" wave function of 
= 1,2) between the ith quasi-a 
3.4 fm, and did not impose any 

















21m 



Fig. 1. The density distributions of mean- field- type 
structure(/3-constraint), "a<- 16 0" structure(d-constraint) 
and "2a- 12 C" structured-constraint), from left to 
right: x are the centroids of the wave packets. 



the right figure, we constrained the distance di 
cluster and the quasi- 12 C cluster to be d\ = d<i 
constraint on the distance between two quasi-a clusters. We see that the nucleons 
are obviously clustering. 

In the case of "a- 16 0" wave functions, when the constrained distance is long, 
two quasi-clusters go nearly to their ground states. This is comparable to the Brink 
model. When the constrained distance is short, nucleon wave packets are deformed 
and wave function is similar to a mean-field-type state. This is different from the 
Brink model. 

In the case of "2a- 12 C" wave 
functions, we constrained the two 
of "a- 12 C" distances, d\ and efo, 
not "a-a" distance. The "a-a" 
distance was optimized by energy 
variation. The 12 C quasi-cluster 
does not need to have a compact 
3a structure, but can have differ- 
ent structure for different values of 
d-constraint parameters d\ and di- 
An important thing to pay atten- 
tion is that nucleons belonging to 
different quasi-clusters can belong 
to the same cluster after energy 
variation. See FigEl which is the 
result of calculation with the con- 



The center of mass of 
a quasi-cluster 1 



a cluster 




O quasi-a cluster 1 
□ quasi-a cluster 2 
x quasi-fe cluster 



The center of mass of 
a quasi-cluster 2 



broken 
s O cluster 



The center of mass of 
12 C quasi-cluster 



Fig. 2. The positions of wave packet centers of nucle- 
ons: The wave function forms "a- (broken 16 O)" 
structure though we constrained two of "a- 12 C" 
distances. 
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straint of d% = 1 fm, c?2 = 2 fm. It can be seen that there are one a cluster and one 
broken 16 O cluster. The reason is that three nucleons of the first a quasi-cluster and 
one nucleon of the second a quasi-cluster make an a cluster because (i-constraint 
does not constrain the position of each nucleon. 

As we mentioned in the introduction, only by the /3-constraint, it is sometimes 
difficult to create the cluster-like wave functions especially for medium-weight nuclei. 
By the d-constraint, we can now create the wave functions rather easy that have 
cluster aspects. 

§4. The Wave Functions Used in the Present GCM Calculation for 20 Ne 

We created 20 Ne wave functions of various types of structures which include 
mean-field-type and a- 16 structure by /3-constraint, and 2a- 12 C and a- 16 struc- 
tures by d-constraint. Only by /3-constraint, we could not create the wave functions 
of "2a- 12 C" structures. The reason is that 2a- 12 C states have higher energy than 
mean-field-type states in small deformation and than a- 16 states in large deforma- 
tion. In the GCM calculation, we combined linearly these 20 Ne wave functions. 

The energy surfaces of the 
wave functions created by 0- 
constraint and by the "a- 16 0" 
distance constraint are shown in 
FigEl In GCM calculation, we 
used x -marked- wave functions 
in the figure. We did not use 
"a- 16 0" wave functions with 
constrained d < 4.0 fm, because 
for d < 4.0 fm nucleon wave 
packets were deformed and the 
wave functions were not differ- 
ent so much from the wave func- 
tions created by /3-constraint. 
The nucleon wave packets of 
"a- 16 0" wave functions with 
constrained d > 4.0 fm are 
spherical and the wave func- 
tions are similar to the Brink 
model wave functions. As for 
the wave functions created by /3-constraint, there is detailed explanation in Ref.2). 

The characteristics of "2a- 12 C" wave functions created by d-constraint are given 
in Table [ij We adopted two types of constraint, d\ = cfe and d\ < d%. We did 
not constrain the distance d^ between two quasi-a clusters in nether cases. As we 
mentioned in H'2.21 in the case of (i-constraint calculation, the resultant distance 
value after energy variation sometimes differs slightly from the aimed distance value. 
For example, in the first three lines of Table QJ the aimed constraint of d\ = d<i is 
slightly violated in the resultant values of d\ and d2, and also the aimed values of 



the distance of a- 16 0[fm] 



— x — d-constraint(a- 16 0) 


— x — x 


x (3-constraint 
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/ 

/ 

/' 

X 

/ 
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/ / 




x - X 





0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 



P 

Fig. 3. The energy surface of wave functions created by 
/3-constraint and that of a- ie O wave functions cre- 
ated by d-constraint: Upper scale is for d-constraint, 
and lower scale is for /3-constraint. Upper and lower 
scales have nothing with each other. 
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Table I. The "2a- 12 C" wave functions used in GCM : di 2 are the distances between the first and 
the second quasi-a cluster and 12 C cluster, d-j, is the distance between two quasi-a clusters, E + 

is the energy of the projected J n = + state in MeV, and Vls is the spin-orbit part of E + in 
MeV. 

di d 2 d 3 E 0+ Vls P 

di = d 2 spherical 3.09 3.15 2.14 -152.4 -10.4 0.36 

wave 3.20 3.24 1.75 -151.7 -11.6 0.37 

packets 3.32 3.33 1.94 -150.8 -11.9 0.40 

deformed 3T3 3T3 L75 -153.7 ^Kl 055 

wave 3.23 3.24 1.82 -152.7 -7.2 0.55 

packets 3.44 3.45 1.96 -150.8 -8.8 0.57 

d L < d 2 deformed 096 L97 L05 -158.2 ^9 043 a- (broken 16 O) 

wave 1.93 2.45 0.70 -158.6 -4.8 0.44 (broken 8 Be)- 12 C 

packets 1.88 2.96 1.74 -158.3 -4.5 0.45 



d\ = d,2 = 3.2, 3.3, and 3.4 fm for the first, second, and third lines, are slightly 
changed in the resultant values of d± and d^- Among the wave functions obtained 
with the constraint of d\ = di, we adopted the wave functions that have two a 
clusters and a 12 C cluster. Only when we constrained d\ = c#2 — 3.2 fm, we got 
such wave functions. Those wave functions are divided into two types of nucleon 
wave packets, spherical and deformed. One of the largest difference between these 
wave functions is spin-orbit energy. The wave functions with spherical nucleon wave 
packets have larger spin-orbit energy than those with deformed nucleon wave packets. 
It is because the 12 C cluster of the spherical wave packets is closer to the ground state 

whose dominant component is the Ops subshell closed structure in the case of the 

2 

Gogny D1S force. But the 12 C cluster of the deformed wave packets has a triangular 
configuration of three prolate deformed a in 12 C cluster which has small spin-orbit 
energy. As the three a in the 12 C cluster are deformed, spin-orbit energy is not zero 
but its magnitude is very small compared with the spin-orbit energy of the ground 
state of an isolated 12 C in the Gogny D1S force which amounts to — 17MeV. As for 
d\ < d-2, we adopted the wave functions that have the lowest intrinsic energy within 
the same condition of constraint. The a- (broken 16 O) wave function has an a cluster 
and a broken 16 cluster by the reason explained in The (broken 8 Be)- 12 C wave 
function has 12 C cluster and eight nucleons not clustering. The last wave function 
obtained with d\ < di has a (2a- 12 C)-like structure, but here also the nucleon wave 
packets are largely deformed and the 12 C cluster has a three deformed-a structure. 

We expected at first that the /3-constraint wave functions in GCM calculation 
would mainly contribute in reproducing K 7 * = O^Oj" bands, the d\ = d2 "2a- 12 C" 
wave functions in reproducing K w = 0% band, and the d\ < cfo "2a- 12 C" wave 
functions in reproducing K n = 0^ band. However, as an actual result, the d\ = di 
wave functions with deformed wave packets and the d\ < di "2a- 12 C" 
wave functions became mainly the components of K w = Of band. 
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§5. Results 



5.1. The energy spectrum 

After the angu- 
lar momentum pro- 
jection, we have 
performed the GCM 
calculation, and ob- 
tained the level scheme 
of 20 Ne(FiggJ). We 
have obtained three 
K n = 0+ bands 
and a K w = 2 + 
band. 

We consider that 
the first, second, 
and third J w = 
+ states of our 
calculations corre- 
spond to the ob- 
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Fig. 4. The excitation energies of the low-lying states of 20 Ne 



served J 77 







and 0j states, respectively. The reason is discussed in the 



next subsection, ^15.21 

5.2. Analyzing of the Results of GCM Calculation 

5.2.1. Radii, deformations, and harmonic oscillator quantum numbers 

The radii, deformation, and harmonic oscillator 
number of the J n = + states are given in Table 
im The harmonic oscillator number operator, N, is 
defined in following way: 



N 



A / ~2 -1 

y ( PL + t mr\ 
^ V 2m 2 * 



/fcj - -A, 



(5-1) 



Table II. Radii, deformations, 
and harmonic oscillator quan- 
tum numbers: R is r.m.s. 
radius in fm, /3 is quadrupole 
deformation parameter, and 
N is harmonic oscillator 
quantum number. 



r 


R 





N 


ot 


3.00 


0.50 


21.2 


ot 


3.22 


0.69 


24.8 


ot 


3.61 


0.87 


31.6 



where A is a mass number, m is the nucleon mass, i p i 
are momentum operators, fj are position operators, 

and is determined so that (N) takes the 

\ / j*-=o+ 

lowest value. 

If wave function is of (s<i) 4 structure, iV is equal 
to 20, and if wave function is of (p)~ 4 (sd) 8 structure, N is equal to 24. The observed 
K w = O2 band is regarded as having a (sd) 4 shell structure, so we consider that 



the second J" 
experiments. 



+ of our calculation does not correspond to the J 71 " = 0t of 
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5.2.2. The components of a- 16 0, 2a- 12 C, and (a- 16 0)U(2a- 12 C) structure 

The amount of the com- 
ponents of a- 16 0, 2a- 12 C and m . . TTT „ r 16 „ „ 12 „ 
v 1fi 19 x Table III. The components of a- O, 2a- C, and 

[a- b O)U(2a- C) structures are (a- 16 0)U(2a- 12 C) structure: E is total energy 

given in Table EH The amount of m MeV, c (a} , c {2a} and c {a}u{2a} are the per- 

the component of {X} structure, centages of the component of q- 16 0, 2a- 12 C and 

C{ X }, is defined as (a- 16 0)U(2a- 12 C) structure, respectively. 



r 


E 




C{2a} 


c {a)U{2a) 


K* 


0+ 


-160.6 


59 


41 


78 


ot 




-152.9 


51 


29 


90 


ot 




-150.6 


69 


14 


92 


ot 


2+ 


-159.6 


58 


31 


72 


ot 




-152.0 


51 


26 


83 


ot 




-150.1 


60 


29 


93 


ot 



($|P X |<£) , V 

where Px is the projection oper- 
ator which projects out the {X} 
component from the $, and Px is 
defined in Appendix. As the basis 
states which span the a- 16 func- 
tional subspace, we adopted the eight wave functions created by the ^-constraint 
on the "a- 16 0" distance d = 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 7.5 and 8.0 fm, as the 
basis states which span the 2a- 12 C, we adopted three wave functions created by the 
d-constraint on the two of "a- 12 C" distance, dip, d\ = c?2 = 3.2, 3.3 and 3.4 fm that 
have two a clusters and a 12 C cluster with spherical nucleon wave packets, and as 
for the combined cluster subspace of (a- 16 0)U(2a- 12 C), we adopted both basis wave 
functions. 

When we compare the second and third J 71 = + of our calculation, we find 
the following points. The component of 2a- 12 C in the second J 71 = + is larger 
than that of third J 77 = + , and the component of a- 16 in the second J n = + 
is smaller than that of third J T = + . This leads that the second J w = + of our 
calculation corresponds to the observed J n = 0^ and the third J n = + of our 
calculation corresponds to the observed J n = 0^ because as explained in ^ the 
observed K 11 = 0^~ states have been considered to have 2a- 12 C structure and the 
observed K n = 0^ states have long been known to have a- 16 structure. 

The percentages of the a- 16 and 2a- 12 C components contained in the calculated 
second J 71 = 2 + state (2^") are similar to those of the calculated second J n = + 
state (Oj), but those of the calculated third 2 + state (2jj~) are slightly different from 
those of the calculated third + state (O3 ), namely the a- 16 component slightly 
decreases from 0^" to 2^" while the 2a- 12 C component slightly increases from 0jj~ to 
2g . Although the 2a- 12 C component of 2% is now comparable to or even slightly 
larger than that of 2^, we classify the calculated 2^~ and 2g" states to the K 71 = 0^ 
and 0^J~ bands, respectively, because the a- 16 component of 2jj~ is larger than that 
of 2^" which, we consider, means that 2jj~ corresponds to the K n = O4 band member. 

We can say that the most importance character of the K n = 0^ band is that it 
has a prominent a- 16 cluster structure and almost all the calculations have never 
failed in reproducing the K n = 0^ band as far as the model space contains the a- 16 
cluster subspace. Therefore we think that our assignment of the calculated 0jj~ and 
2^j~ states to the K n = 0^ band members is rather reliable. 

On the other hand, we think the our present treatment of the 2a- 12 C subspace 
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is not sufficient, because, as we explained in fjfl the GCM basis states which can 
be said to have 2a- 12 C structure are only three states obtained by the d-constraint 
with d\ = 0I2 = 3.2, 3.3, and 3.4 fm. If we include more basis states with 2q- 12 C 
structure, the percentages of the 2a- 12 C component in the calculated Og" and 2^" will 
increase while those in the calculated and 2^~ will decrease. 

5.2.3. The components of the binding energy 

The decomposition of 
the calculated binding en- 
ergy is given in Table IIVI 



J 71 " = O4 state has less spin- 
orbit energy than J 71 " = 0^ 



Table IV. The components of the binding energy: E is total 
energy, T is kinetic energy, V2 is 2-body energy, Vls 
is spin-orbit energy, Vc is Coulomb energy, and V g is 
density dependent energy in MeV 



and O3 states 







It is because 
state has more 



a- 16 component, of which 
spin-orbit energy is nearly 
equal to zero, than J w = 0^ 
and states. 



r 


E 


T 


V 2 


Vls 


V c 


v s 


of 


-160.6 


283.5 


-892.0 


-5.7 


19.1 


434.5 


ot 


-152.9 


280.2 


-866.2 


-5.0 


18.6 


419.5 


oi 


-150.6 


273.3 


-847.9 


-3.2 


17.8 


409.3 



§6. Summary 



We have proposed a new constraint of clustering for the AMD calculation, which 
we call d- constraint. By the d-constraint, we can create rather easily the wave 
functions that have cluster aspect. It is very helpful and important for the AMD 
study because in AMD approach we sometimes find difficulty in getting the wave 
functions with some specific cluster structure. 

We have applied the d-constraint to 20 Ne nucleus. We could create rather easily 
2a- 12 C wave functions which were not easy to get in the AMD approach. Using 
these wave functions, we could reproduce K 77 = 0jj~ band which contains the 2a- 12 C 
structure as an important component. It should be noted that even the existence 
of the K n = 0jj~ band has not been shown before not only by AMD but also by the 
Hartree-Fock approach. 

Now we can create rather easily the cluster-like wave functions for medium- 
weight nuclei. We expect that our new constraint is powerful for the study of the 
coexistence of the cluster structure and mean-field structure in sd and p/-shell nuclei 
and the study of the relation between superdeformation and cluster structure. 
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Appendix 



Suppose the states \4>i) span the functional space {X}. Their overlap matrix By- 
is defined as 

B, 



'13 — \Yt\ YJ/ 



The orthonormalized basis wave functions of the {X}, 
combination of the \4>i), 



I Pa 



(Al) 

, are given by the linear 
(A2) 



where the p a and Cj a are the eigenvalues and eigenvectors of the overlap matrix Bij , 
respectively, 

3 

The eigenvectors Ci a form a complete orthonormalized set, 



c ia c ja 



(A-4) 



It is easy to show that 



are orthonormal. 



Using 



a), Px is defined as, 



a / \ 



(A5) 
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